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We report on the experimental observation of thermoelectric currents in superconductor- 
ferromagnet tunnel junctions in high magnetic fields. The thermoelectric signals are due to a 
spin-dependent lifting of particle-hole symmetry, and are found to be in excellent agreement with 
recent theoretical predictions. The maximum Seebeck coefficient inferred from the data is about 
— 100 pV/K, much larger than commonly found in metallic structures. Our results directly prove 
the coupling of spin and heat transport in high-held superconductors. 


The coupling of electronic and thermal properties in 
nanoscale superconductor hybrid structures is an active 
field of current research due to applications as thermome¬ 
ters, microrefrigerators and particle detectors [l|. Re¬ 
cently, large thermoelectric effects were predicted to oc¬ 
cur in superconductor-ferromagnet hybrid structures 
1^. Thermoelectric effects in metals are driven by bro¬ 
ken symmetry between electron and hole carriers and 
are usually quite small, of the order of a few p.V/K at 
room temperature and vanishing at low temperature. 
While superconductors obey overall electron-hole sym¬ 
metry, spin-splitting of the quasiparticle density of states 
breaks electron-hole symmetry for each spin band, and 
in conjunction with the spin-dependent conductance of 
a superconductor-ferromagnet tunnel junction leads to 
thermoelectric effects Q . The characteristic energy scale, 
the energy gap of the superconductor, is very small com¬ 
pared to the Fermi energy, and consequently, the thermo¬ 
electric effects are predicted to be large. Apart from be¬ 
ing large, which makes them promising for applications, 
these thermoelectric effects are also interesting from a 
fundamental point of view. They lead to a coupling of 
spin and heat currents, which forms the basis of spin 
caloritronics . Thermally driven spin currents have re¬ 
cently been proposed Q as an explanation for the exper¬ 
imentally observed long-range quasiparticle spin trans¬ 
port in superconductors 0, 0. Here, we report on the 
experimental observation of thermoelectric currents in 
high-held superconductor-ferromagnet tunnel junctions. 

The current through a tunnel junction in the presence 
of a voltage V and a temperature difference ST across 
the junction can be conveniently described in the linear 
regime by 

I = gV + VY 

where g is the conductance, T is the average temperature, 
and r] describes the thermoelectric current. It is related 
to the Seebeck coefficient S = —VjdT measured in an 
open circuit by ?; = SgT. Measuring the thermoelectric 
coefficient g is the main purpose of this work. 

Our samples were prepared by e-beam lithography and 



FIG. 1. (color online) (a) False-color scanning electron mi¬ 
croscopy image of a region of sample FISl with measure¬ 
ment configuration for the thermoelectric measurements, (b) 
Sketch of the measurement conhguration for conductance and 
heater calibration measurements. 

shadow evaporation technique. First a thin supercon¬ 
ducting (S) aluminum strip of thickness ^ai ~ 20 nm was 
evaporated, which then was oxidized in situ to form a 
thin insulating (I) tunnel barrier. After that ferromag¬ 
netic (F) iron (tpe ~ 15 — 20 nm) and normal-metal (N) 
copper wires (tcu ~ 50 nm) were overlaid from different 
angles. In this manuscript we present and compare re¬ 
sults from three samples, two samples with ferromagnet- 
superconductor junctions (FISl and FIS2), and a non¬ 
magnetic reference sample (NIS). Figure [Ha) shows a 
false-color scanning electron microscopy image of a re¬ 
gion of sample FISl. The central part of the structure 
is a 6-probe junction (a) consisting of a superconductor- 
ferromagnet tunnel contact overlaid with an additional 
copper wire. This design allows us to measure the cur¬ 
rent through the tunnel junction while simultaneously 
passing a heater current along the iron wire to create a 
temperature difference across the junction. In addition, 
there is a 4-probe normal-metal junction (/3), located at 
a distance d w 1.5 /im from a. This junction is used for 
control experiments. Sample FIS2 has a slightly different 
layout (not shown) which avoids passing the heater cur¬ 
rent across the junction. As it turns out this had no effect 
on the results. The third sample (NIS) has the same lay¬ 
out as FISl, with the iron wire replaced by a copper wire 
of a similar thickness. The measurement schemes are de¬ 
scribed below together with the results. Measurements 
were performed in a dilution refrigerator down to a base 
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FIG. 2. (color online) (a) Differential conductance g of junc¬ 
tion a of sample FISl as a function of bias V for different 
magnetic fields B. (b) Zero-bias conductance g as function of 
magnetic field B for different base temperatures Tq. (c) Differ¬ 
ential conductance g at fixed base temperature and magnetic 
field for different heater currents /heat- (d) Temperature Tp 
of the ferromagnet as a function of /heat for different Tq. 


temperature of Tq = 50 mK, with an in-plane magnetic 
field B applied parallel to the iron wire as indicated in 
Fig.HKa). We denote the temperature of the refrigerator 
by To to distinguish it from the electronic temperature 
of the sample throughout this manuscript. 

Before we describe the thermoelectric measurements, 
we first characterize the sample and calibrate the heater 
current. In Figs. [2Ka) and (b) we show the differential 
conductance g of junction a of sample FISl, measured 
by standard low-frequency ac technique in the configura¬ 
tion shown in Fig. [IJb), without heater current. Figure 
Ha) shows 5 as a function of bias voltage V for different 
magnetic fields B. In high fields the Zeeman-splitting of 
the density of states is visible together with a broadening 
of the curves due to orbital pair breaking. To fit our data, 
we model the current using eqn. (2a) of Ref. d. We as¬ 
sume that the superconductor is at temperature Ts = T, 
and the ferromagnet is at temperature Tp = T+6T, with 
voltage V applied to the ferromagnet. This yields 
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where Gt is the normal-state junction conductance, P 
is the spin polarization of the junction conductance, 
e = — |e| is the charge of the electron, and /o is the Fermi 


function. Equation (ED includes conductance (V ^ 0, 
5T = 0) and thermoelectric currents {V = 0, ST ^ 0) 
on an equal footing. The density of states factors are 
Nq = (7V+ -|- 7V_)/2 and = iV+ — iV_, where N± 
are the densities of states for the two spin projections in 
the superconductor. and fV_ were obtained from the 
standard model of high-held tunneling 00 , including 
the effect of Zeeman splitting, orbital pair breaking and 
spin-orbit scattering. The factor is odd in energy, 
and gives rise to the observed thermoelectric effect. It 
is non-zero only in the presence of a Zeeman splitting of 
the density of states in combination with a spin-polarized 
junction conductance. Fits of our data based on eqn. ED 
are shown as lines in Fig. Ha)- From these fits we extract 
the spin polarization P « 0.08. Due to the finite spin po¬ 
larization P, there is a small but visible asymmetry of the 
conductance at high fields [HI , with higher conductance 
for electrons (negative bias) than holes (positive bias). 

In Fig.Hb) we show g at zero bias as a function of the 
magnetic field B for three base temperatures Tq, together 
with fits to our model. For these fits, the self-consistent 
pair potential A, orbital pair-breaking strength Oorb and 
effective Zeeman splitting were calculated as a function 
of the magnetic field using Ref. |I3, including the effect 
of Fermi liquid renormalization on the effective Zeeman 
splitting. The latter, with the Fermi liquid parameter 
G° = 0.3 taken from literature [0, was found to im¬ 
prove the fits close to the critical field. The fits of the 
conductance are used to establish the necessary sample 
parameters for modeling the thermoelectric results shown 
below. 

We now turn to the heater calibration. In Fig. He) 
we show the differential conductance 5 as a function of 
bias V at fixed temperature Tq = 100 mK and mag¬ 
netic field P = 1 T for different dc heater currents /heat 
(see Fig.Hb) for the schematics). Upon increasing /heat, 
the conductance visibly broadens due to increasing tem¬ 
perature Tp of the ferromagnet (the thermal broadening 
depends on Tp only and is independent of Ts). From fits 
of the data we extract Tp as a function of /heat, shown 
in Fig. Hd) for three different base temperatures Tq. 

To describe the heating of the junction as a func¬ 
tion of current, we assume quasiequilibrium with negligi¬ 
ble electron-phonon scattering, as appropriate for meso¬ 
scopic metal wires of about 10 /rm length at sub-Kelvin 
temperatures [l|. In this case, the temperature at the 
junction (in the middle of the heater wire) is given by 


Tp 
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where //heat is the resistance of the heater wire, Lq = 
TT^fcg/Se^ is the Lorenz number, and T is the electronic 
base temperature in the absence of heating. In Fig. Hd) 
we show fits of our data using eqn. (H- As can be seen, 
the fits show a good agreement with the data, and we 
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FIG. 3. (color online) (a) Differential conductance g of con¬ 
tact P of sample FIS2 at fixed magnetic field B and temper¬ 
ature To for different heater currents /heat- (b) pair potential 
A as a function of /heat obtained by fitting the differential 
conductance shown in (a), (c) increase 5Ts of the tempera¬ 
ture of the superconductor as a function of /heat extracted by 
inverting A{Ts), together with STf. 


use the fits for temperature calibration for the thermo¬ 
electric measurements described below. The electronic 
base temperature T is slightly increased over the cryo¬ 
stat base temperature Tq at low temperatures, proba¬ 
bly due to incomplete filtering of thermal noise from 
higher temperature stages of the cryostat. The fits yield 
Rhea.t ~ 230 n, smaller than the 2-probe resistance of the 
iron wire, //pe = 1350 fl. We attribute this to the fact 
that the thick copper wire acts as a cooling fin. 

Superconductors are poor heat conductors at low tem¬ 
perature, and can therefore be easily heated. To check 
the impact of the heater current on the temperature Ts 
of the superconductor we performed additional control 
experiments on junction j3. The quasiparticle energy re¬ 
laxation length is typically a few 10 /rm in aluminum at 
low temperatures 1^, and we assume Ts to be nearly 
the same at contacts a and /3. An increase of Ts affects 
the differential conductance g only indirectly by a small 
reduction of the pair potential A. In Fig.[3]we show data 
of sample FIS2, for which we made the most detailed 
measurements of Ts. Figure [3Ka) shows the differential 
conductance g for different heater currents /heat at fixed 
magnetic field // = 1.3 T and temperature Tq = 100 mK. 
As can be seen, there is almost no change at small cur¬ 
rents, and at larger currents, the gap is slightly reduced. 
The pair potential A obtained by fitting the data is plot¬ 
ted in Fig. IlKb) as a function of the heater current. The 
change of A is of the order of a few /xeV for the largest 
applied current, with considerable scatter on this small 
scale. To reduce scatter, we averaged data for two adja¬ 
cent points and then inverted the self-consistent relation 
A(T -I- STs, B) to obtain dTs (assuming dTs = 0 for the 
first point). In Fig. [SKc) the resulting 5Ts is shown as 
a function of /heat, together with dTp. For thermal bias 
(5Tp < 100 mK used in the thermoelectric experiments 
(indicated by the dashed line), we find 5T^ < 20 mK, 
much smaller than JTp. We conclude that most of the 



FIG. 4. (color online) (a) Thermoelectric current Ah as 
a function of magnetic field B at base temperature To = 
100 mK, together with theoretical calculations (see text), (b) 
thermoelectric coefficient g of sample FISl for different tem¬ 
peratures To. (c) thermoelectric coefficient ?7 at To = 100 mK 
for all three samples, (d) Seebeck coefficient S' as a function 
of B calculated as described in the text. 


thermal bias actually drops across the junction. 

Now we will focus on the thermoelectric measurements. 
For these measurements, we applied an ac heater cur¬ 
rent, as shown schematically in Fig. Ha). Since the 
heating power is proportional to /peat; thermal ex¬ 
citation appears on the second harmonic of the signal, 
and we detected the second harmonic of the resulting 
current Ah flowing into the superconductor. The heater 
current amplitude was typically chosen to yield a peak- 
to-peak thermal amplitude STp = 100 mK. In Fig. |D(a) 
we show the peak-to-peak amplitude of the thermoelec¬ 
tric current Ah as function of the magnetic field B at 
base temperature Tq = 100 mK for sample FISl. As 
expected, there is no signal at zero field, be., in the ab¬ 
sence of Zeeman splitting. With increasing field a neg¬ 
ative signal develops above i? = 0.5 T and reaches a 
broad maximum around i? = 1 T, then going back to 
zero as B approaches the critical field. The thermoelec¬ 
tric current is negative, as expected for our junction with 
higher conductance for electrons than holes, as seen in 
Fig. [2Ka). We also show theoretical calculations using 
eqn. ([2]) with the same parameters as for the fits of the 
conductance shown in Fig. [2jb). For the amplitude 5T 
we chose two assumptions: The dashed line is calculated 
using 5T = 5Tp = 100 mK obtained from the calibration 
fit to eqn. whereas the dash-dotted line is calculated 
using ST = STp — STs ~ 80 mK. The data lie in between 
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both curves, and the best fit (solid line) is obtained for 
6T = 86 mK. In Fig. HKb) we show the thermoelectric 
coefficient rj = I^i^T/dT for sample FISl as a function of 
the magnetic field B for different base temperatures Tq. 
Here T = T + ST(2, and we used ST from the best fits. 
For increasing temperatures we see that both the mag¬ 
nitude and broadening of rj increase, as expected from 
theory Q. Figure SKc) shows 77 as a function of the nor¬ 
malized field B/Bc for all samples at Tq = 100 mK. The 
data are similar for both ferromagnetic samples FISl and 
FIS 2. There is no signal for the non-magnetic reference 


directly prove the coupling of spin and heat transport in 
high-held superconductors. 


Present addresss: Institute of Technical Physics, Karl¬ 
sruhe Institute of Technology, Karlsruhe, Germany 
e-mail address: detlef.beckmann@kit.edu 
F. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin. 


sam ple NIS. directly proving the spin-dependent origin 


[ 1 ] 

and J. P. Pekola, [Rev. Mod. Phys. 78, 217 (2006) 

[21 P. Machon. M. Esch r ig, and W. Belzig, jPhys^RewLett 


110, 047002 (2013); New J. Phys. 16, 073002 (2014) 


of the effect. From the data and Hts in Figs. Illa)-(c), we 
conclude that the data are in excellent qualitative and 

guaiitilHlive agieeiiieiit with the llieuietical ptedicliuiis 

- 


[3] A. Ozaeta. P. Virtanen. F. Bergeret . and T. Heikkila, 


Phys. Rev. Lett. 112, 057001 (2014) 


[41 M. S. Kalenk o y and A. D. Zaikin, Pl^s^^^ew^ 


90 


134502 (2014); Phys. Rev. B 91, 064504 (2015) 


Ii j Fig.Hfd) we hnally show the Seebeck coefficient S = Mater. 11, 391 (2012) 


|5| G. E. W. Bauer, E. S aitoh, and B. J. van Wees, Nat 


r]/{gT) inferred from our measurements as a function of 
the applied magnetic held B for sample FISl. The data 
and hts of S were both calculated from the data and hts 
of 77 and g shown in Figs. SKb) andUKb). As can be seen, 
S shows a similar held dependence as 77 itself, however 
with a faster increase and maximum at smaller helds. 


[6J M. Silaev, P. Virtanen, F. S. Bergeret, and T. T. 


and 


Heikkila, Phys. Rev. Lett. 114, 167002 (2015). 

[7] F. Hiibler, M. J. Wolf, B. Reckmaun. 

H. V. Lohneysen, Phys. Rev. Lett. 109, 207001 (2012); 
M. J. Wolf, F. Hiibler, S. Kolenda, H. v. Lohneys en. and 


The Seebeck coefficient is larger at lower temperature 


D. Beckmann 
Wolf, G. SiirEers. G. Fisch er, 


Phys. Rev. B 87, 024517 (2013); M. J. 


Rev. B 90, 144509 (2014) 


and D. Beckmann, Phys 


due to the freeze-out of quasiparticle conductance. As a 
result, the scatter is also amplihed, and the data for Tq = 
100 mK are unreliable and omitted from the plot. For 


Nat. Phvs. 9. 84 120131 


]9] K. Maki, 

Prog. Theor. 

Phys. 32, 29 (1964) 


and M. Aprili, 


To =|: 250 mK, the maximum value Is 5” ~ —100 /iV/K. 


[101 R. Meservev, P. M. Tedrow, and R. G. Bruno, Phys. Rev 


B 11, 4224 (1975) 


I n conclusion, we have measured thermoelectric cur- [111 P. M. T edrow and R. Meservey, Phys. R ev. Lett. 26, 192 


higi 


5 in siioerconductor-ferromaffnet tunnel iunctions at 

(1971); R. 

magnetic fields. The results are in excellent agree- 

173 (1994) 


R. Meservey and P. M. Tedrow, Phys. Rep. 238 


ment with recent theoretical predictions. The Seebeck 
coefficients inferred from the data can be as large as 
100 77 V/K, much larger than usually found in metal 
structures, despite the fact that the spin polarization of 
our tunnel junctions is relatively small. With larger spin 
pol ariza,tion, e.g., u.sing .spin-filter junctions [T^. Seebeck 


[12J J. A. X. A 


exander, T. P. Orlando, D. Rainer, and P. M. 


Tedrow, Phys. Rev. B 31, 5811 (1985) 

[13] P. M. Tedrow. J. T. Kucera. D. Rainer , and T. P. Or¬ 
lando, Phys. Rev. Lett. 52, 1637 (1984)| . 

[14] K. Y. Arutvunov. H.-P. Aurane va. and A. S. Vasenko, 


Phys. Rev. B 83, 104509 (2011) 


[151 X. Hao. J. S. M oodera. and R. Meservey, Phys. Rev. B 


coe fficients exceeding 1 mV/K seem feasible. Our results 


42, 8235 (1990) 





































































